clear all
set more off
capture log close
log using "$output/figureA2.log", replace

version 16	

use "$data/ctctaxsim_mar1kid1.dta", clear
append using "$data/ctctaxsim_mar1kid2.dta"
append using "$data/ctctaxsim_mar1kid3.dta"

*** 2017 law
twoway line v25 earnings if year==2017 & depx==1 || line totctc earnings if year==2017 & depx==1, graphregion(color(white)) bgcolor(white) text(2.600 18.000 "EITC" , j(left) size(medsmall)) text(.200 25.000 "CTC" , j(left) size(medsmall)) ylabel(0(5.000)10.000, axis(1)) leg(off) ytitle("Tax Credits ($1000s)") xtitle("Earnings ($1000s)" "{bf:Panel A: 2017 Policy}") text(9.500 80.000 "2017 Policy" , box fc(white) margin(small) j(left) size(medsmall)) || line taxesowed earnings if year==2017 & depx==1, yaxis(2) ylabel(0 "0" -5.000 "5" -10.000 "10", axis(2)) lp(-) text(3.500 75.000 "Income Tax" "(Right Axis)" , j(left) size(medsmall)) ytitle("Income Taxes Owed ($1000s)", axis(2))
graph save "$output/ctceitc1_mar1kid1.gph", replace

*** 2018 law
twoway line v25 earnings if year==2018 & depx==1 || line totctc earnings if year==2018 & depx==1, graphregion(color(white)) bgcolor(white) text(2.800 18.000 "EITC" , j(left) size(medsmall)) text(.500 22.000 "CTC" , j(left) size(medsmall)) ylabel(0(5.000)10.000) leg(off) ytitle("Tax Credits ($1000s)") xtitle("Earnings ($1000s)" "{bf:Panel B: 2020 Policy}") text(9.500 80.000 "2020 Policy" , box fc(white) margin(small) j(left) size(medsmall)) || line taxesowed earnings if year==2018 & depx==1, yaxis(2) ylabel(0 "0" -5.000 "5" -10.000 "10", axis(2)) lp(-) text(4.500 75.000 "Income Tax" "(Right Axis)", j(left) size(medsmall)) ytitle("Income Taxes Owed ($1000s)", axis(2))
graph save "$output/ctceitc2_mar1kid1.gph", replace

*** 2021 law
twoway line v25 earnings if year==2018 & depx==1 || line ctc2021 earnings if year==2018 & depx==1, graphregion(color(white)) bgcolor(white) text(4.000 20.000 "EITC" , j(left) size(medsmall)) text(3.500 35.000 "CTC" , j(left) size(medsmall)) ylabel(0(5.000)10.000) leg(off) ytitle("Tax Credits ($1000s)") xtitle("Earnings ($1000s)" "{bf:Panel C: 2021 Policy}") text(9.500 80.000 "2021 Policy" , box fc(white) margin(small) j(left) size(medsmall)) || line taxesowed earnings if year==2018 & depx==1, yaxis(2) ylabel(0 "0" -5.000 "5" -10.000 "10", axis(2)) lp(-) text(5.00 75.000 "Income Tax" "(Right Axis)" , j(left) size(medsmall)) ytitle("Income Taxes Owed ($1000s)", axis(2))
graph save "$output/ctceitc3_mar1kid1.gph", replace

*** THIS COMPARES TAX CREDITS IN 2017, 2018, 2021
twoway line negfiitax earnings if year==2017 & depx==1, lp(-) pstyle(p1) || line negfiitax earnings if year==2018 & depx==1, lp(_) pstyle(p2) || line negfedtax2021 earnings if year==2018 & depx==1,  pstyle(p1) graphregion(color(white)) bgcolor(white) ylabel(-10.000 -5 0 5 10.000) leg(order(1 "2017" 2 "2020" 3 "2021") c(1) pos(7) ring(0)) ytitle("Net Tax Credits ($1000s)") xtitle("Earnings ($1000s)" "{bf:Panel D: Net Effect from CTC, EITC, Income Tax}") text(5.00 75.000 "Net Tax Credits" "from EITC, CTC," "Income Taxes", box fc(white) margin(small) j(left) size(medsmall))
graph save "$output/ctceitc4_mar1kid1.gph", replace


graph combine "$output/ctceitc1_mar1kid1.gph" "$output/ctceitc2_mar1kid1.gph" "$output/ctceitc3_mar1kid1.gph" "$output/ctceitc4_mar1kid1.gph" ,graphregion(color(white)) 
graph export "$output/figureA2.png", replace

erase "$output/ctceitc1_mar1kid1.gph"
erase "$output/ctceitc2_mar1kid1.gph"
erase "$output/ctceitc3_mar1kid1.gph"
erase "$output/ctceitc4_mar1kid1.gph"

log close